arXiv: 1501.06617vl [hep-lat] 26 Jan 2015 


PROCEEDINGS 

OF SCIENCE 



Tuning of the strange quark mass with optimai 
reweighting 


Bjorn Leder"^, Jacob Finkenrath 

Department of Physics, Bergische Universitdt Wuppertal 

Gaussstr. 20, 42119 Wuppertal, Germany 

E-mail: leder@physik . uni-wuppertal. de 


Quark mass reweighting can be used to tune the mass of dynamical quarks. The basic idea is 
to use gauge field ensembles generated at some bare mass parameters to evaluate observables at 
different bare sea quark masses. This involves the computation of so called reweighing factors 
which are given as ratios of fermion determinants. In the case of simulations including the strange 
quark, reweighting can be used to improve the approach towards physical quark masses. Optimal 
reweighting strategies combine a change of the strange quark mass with a change of the light quark 
masses in order to minimize the fluctuations of the reweighting factor. We present numerical test 
of such strategies for recent CLS2 simulations and a software package for mass reweighting based 
on openQCD. 
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1. Why mass reweighting? 


Lattice QCD simulations proceed to generate ensembles at small lattice spacings and light 
quark masses close to or at their physical values. This is possible because of the algorithmic im¬ 
provements of the last decade. Nevertheless, the cost to generate independent gauge conhgurations 
grows, particularly so because of the growing autocorrelation times when the lattice spacing is low¬ 
ered [1]. The problem is further aggravated by the inclusion of a sea strange quark, which leads to 
an enlarged mass tuning problem (see for example [2]). 

In such a situation it might be more cost efficient to change the masses of a given ensemble 
using reweighting instead of generating a new one. This is because although reweighting is ex¬ 
pensive it only has to be done for independent gauge conhgurations and therefore its cost does not 
increase with the autocorrelation times. Also the naively expected linear scaling with the volume 
might be much milder in reality [3, 4]. 

Reweighting simply means changing the probability Pa{U) (“weight”) of each gauge conhgu- 
ration U in the expectation value 

{^)a = Y J^[U]Pa{U)^{U), Pa{U)=e-^^‘^l^’^^p_dct{D{U)+mi + ir5Pi). ( 1 . 1 ) 


Here we assume Hf quark havors with (possibly degenerate) bare masses nt; and twisted masses 
Pi, the gauge action Sg only depends on j3 and normalization /D[f/]P(f/)/Z = 1. The subscript 
on the expectation value and the probability distribution specihes the bare parameter set a = 
{j3,mi,ni2,... ,mn^.,pi,... ,Pnj,}- Now an observable at bare parameter set h = {p' 
ju(,...,is obtained from 


{^)b = 




( 1 . 2 ) 


If one restricts oneself to changes in the masses (hence p' = p), as we do here, the reweighting 
factor W can be written as 


1 

W = - 

det(A) ’ 


«/ 

^=n 


r=l 


D + mj + ijsPi 
D + nf- + ijsp'i ' 


(1.3) 


The next section specihes the stochastic estimator used for the determinant in (1.3). In section 
3 follows a summary of the properties of three important cases of mass reweighting. From these 
properties an optimal strategy for the tuning of the strange quark mass is derived in section 4. 
Section 5 briehy describes a publicly available implementation of this strategy and in section 6 
numerical tests and cost estimates are presented. 


2. Fluctuations 


The numerical estimation of reweighting factors of the form (1.3) uses the Monte-Carlo eval¬ 
uation of an integral representation of the determinant of a complex matrix [5, 6] 
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( 2 . 1 ) 
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with N-q Gaussian distributed random vectors T](') and -^W for —)■ oo if A(A+A^) > 0 
[5, 6]. Let us assume A can be written as A = / + eB with e| |B| | ^ 1. Then the stochastic error of 
this estimator may be expanded 

52(A) = var,(l%J/(iV,|lT|2) = + (2.2) 

rs^ 

A similar expansion is obtained for the fluctuations of the reweighting factor in the expectation 
value ( 1 . 1 ) 

ol = varc/(lT)/ (IT)^ = Ehwu{Tv{B)) + 0{e^). (2.3) 

It has been noticed that the stochastic estimator (2.1) may be plagued by large fluctuations 
and/or long tails in the distribution. The reason is that its variance may not be defined even if the 
mean is, i.e., the variance is only defined for A (A +A^) > 1 [7]. There have been several attempts 
to control the variance [ 8 , 9, 10]. Here we use a factorization of the determinant in N factors with 
e' oc e/N [5, 6 ]. Since the factorization has been described thoroughly elsewhere and for better 
readability we skip the details here and stick to the formulas with N = 1. The generalization is 
straightforward. 


3. Mass reweighting factors 


3.1 One flavor 


The simplest case of mass reweighting is the change of the mass of one quark flavor of mass 
nis to mass nis + Am. Using the shorthand D,n = D + m, the corresponding reweighting factor is 
(Bs = 0 ) 


Wif{ms,Am) = 


det (Z),jjj+Am) 


Aif = /-AmD ' 


det(D,„J det(A,f)’ ‘ 

and the stochastic error and the fluctuations scale (asymptotically) with Arn^. 


(3.1) 


3.2 Two flavors 


Now we consider the simultaneous change of the mass of two quark flavors s and 1. Without 
loss of generality we assume ni/ < nis and a change mj —)■ ni/ — yAni while nis —^ ms + Am. The 
corresponding reweighting factor is (/r^ = Bi = 0) 


W2p{mi,ms,Am) = 


det{D,n, —yAmDirts+Am 


det(A 2 f) ’ 


A 2 f = / + Am 


YAm + YD,„^-D„y 
^nts +Am^mi—yAm 


(3.2) 


and the stochastic error and the fluctuations scale (asymptotically) with Ain^. However, if m/ = 
nis = m and 7 = 1 one finds A 2 f = / + Am^(D^ — Am^)^^ and a scaling oc Am^ (see the isospin 
reweighting in [ 6 , 4]). In this case the noise and the fluctuations of decreasing the mass of one 
quark are compensated by increasing the mass of the second one. For 7=0 the reweighting factor 
reduces to the one flavor case: {mi,ms,Ani) = Wif{nis,Am). Therefore, in the general case 

nil < ms, an optimal 0 < 7* < 1 may be found that minimizes the fluctuations of . With light 
quarks around the strange quark mass an optimal value of 7 ^ « 0.82 was found [5]. 

A generalization of eq. (3.2) for finite jXs and/or fXi is straightforward. Note that, in general, in 
this case the reweighing factor is complex. 
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(i) 

(ii) 

(iii) 

(iv) 

up 

(m/,0) 

(m/,/l) 

—)• (m/ — yAm,/^) 


—)• (m/ — yAm,0) 

down 

strange 

(m/,0) 

(m^,0) 

(«^/,-At) 

—>• (m^ + Am,0) 

—)• (m/ — yAm, —/i) 
—>• (m^ + 2Am,0) 

—)• (m/ — yAm,0) 

type 


tm (sec. 3.3) 

2f (sec. 3.2) 

2f (sec. 3.2) 

tm (sec. 3.3) 


Table 1: Optimal reweighting strategy for strange quark mass reweighting. 


3.3 Twisted mass reweighting 


Twisted mass reweighting was introduced [3] and implemented [11] in order to stabilize the 
HMC at small quark masses. In the context of mass reweighting it can be used to ensure the 
convergence of the stochastic estimator (2.1) [5]. For a doublet of quarks of mass m\ = m 2 = m 
and /ii = jU = —/f2 the corresponding reweighting factor is 




det{D^Di + lx'^) _ 1 ^ 

det(D^D,^ + At2) det(Atn>)’ "" D,„Dj, + At'2’ 


(3.3) 


where we used (D,„ + /ys/i)^ = Y 5 {D,„ — /75/i)75. The stochastic error and the fluctuations scale 
(asymptotically) with (/i^ — 


4. Optimal strange quark mass reweighting 

As pointed out in the introduction, inany = 2 + l simulation with up, down, strange quark 
masses {ni/, m;, m^} it might be necessary to tune the strange quark mass of an ensemble in order 
to improve the approach to the physical point. We have seen in section 3 that it is advantageous to 
combine the change of the mass of one quark with an opposite change of another one. Furthermore, 
if light quarks are reweighted, a finite twisted mass serves as a safeguard against zero crossings of 
small eigenvalues. Taken together, these lessons lead us to propose the following reweighting 
strategy: 

(i) the light quarks are reweighted to a finite /i, 

(ii) the up and the strange quark are reweighted together in opposite directions, 

(iii) the down and the strange quark are reweighted together by the same amount, 

(iv) the light quarks are reweighted to zero twisted mass. 

The four steps are summarized in Table 1 . In the stochastic estimation the twisted mass reweighting 
in step (i) and (iv) can be combined so that the total reweighting factor can be written as a 
product of three factors 

W, = W(ii)W(iii)W(i+iv) ^ (W(ii))^„ (W(iii)),v, (W(Wv));v,, (4.1) 

and each of these factors is estimated according to (2.1). The total change in the masses is 
{mi, mi,ms} ^ {mi — yAm, m/ — yAm, + 2Am} . 
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If the ensemble is generated at finite jj. in the light quark sector, as for example in [2], step (i) 
can be omitted. Depending on which kind of twisted mass reweighting is used in the production of 
the ensemble, the last factor is then given by 


W(i+iv) = < 


W^tm(m/, 0, /I) Wtm {mi - Y^m, /i, 0) 

Wiyn{mi-YAm,iJ.,0) 

Wtm {mi , H , \/2jU) Wtm {mi - J^m, jj ,, 0) 


none 

type I in [11] . 
type II in [11] 


(4.2) 


The reweighting factors for step (ii) and (iii) are explicitly given by 

W(ii) = {mi , m,, Am), with ft/ = ft (4.3) 

W(iii) = {mi , ms + Am, Am), with ft/ = —ft, (4.4) 

and 7 *. As mentioned before for ft/ / 0 the reweighting factor is complex. But since the 
product W(ii)W(iii) is manifestly real, the phases of the two factors cancel. 


5. Public code for mass reweighting 

The numerical results presented in the next section have been obtained with the mrw-package 
[12], an extended version of the openQCD package [13]. Keeping the structure of openQCD the 
extension is implemented as a module (mrw). It provides a main program that reads in open¬ 
QCD style input files, documentation and sample input files. The mrw-package is publicly avail¬ 
able at https : //github. com/b joern-leder/mrw. In detail it adds the following features to 
openQCD: 

• one flavor mass and twisted-mass reweighting [5, 6 ] 

• interpolation (factorization) for twisted mass reweighting type I and II 

• factorization with non-equidistant interpolations [5, 6 ] 

• isospin mass reweighting [5, 6 ] 

• strange quark mass reweighting of section 4 

• several check routines for all parts of the new module 


6. Numerical test and cost 

The optimal strange quark mass reweighting proposed in section 4 has been tested on two 
jif = 2+ I ensembles from the effort described in [2], see Table 2. While the first one is off the 
quark mass trajectory described in section 2.3 of [2], the second one sits on it.* The lattice spacing 
is approximated to a = 0.086 fm [2] and the lattice size is 64 x 32^. 

'This trajectory is defined by X!/™/ = const. Note that for 77^ 1 the optimal strange quark reweighting takes the 
ensemble off, whereas for 7 = 1 it would move the ensemble along this trajectory. 
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ID in [2] 

ffiji [MeV] 

mK [MeV] 

Am, [MeV] 

7 


- 

330 (380) 

450(430) 

-12 

0.80 

0.0 

BIOS 

280 ( 200 ) 

460 (480) 

12 

0.80 

0.001 


Table 2: Ensembles used in the numerical test. In parentheses are the meson masses after reweighting. 
AiHi is the difference of the renormalized strange quark masses before and after the reweighting. The first 
ensemble was part of an exploratory study and is not listed in [2]. 



0 50 100 150 200 250 300 350 

MDU 

Figure 1: Strange quark reweighting factor Ws for ensemble B105. 


In Figure 1 the reweighting factor (eq. (4.1)) is plotted for the ensemble B105. The phase 
is also plotted and is compatible with zero as expected. The fluctuations of the reweighting are 
sizable, but the change in the meson masses is large: 5% (30%) for the kaon (pion). By using a 
smaller value for / this change can be made more balanced. Assuming a scaling of the fluctuations 
oc Am^V‘‘/rni the feasibility of such a reweighting at different parameters can be projected using 
the following formula 


/ Am, y ydOMeVy 
Vl2MeV J V ^ ) (3.3fm)4‘? ’ 


^ = 1/4...3/4. 


where Am, is the difference of the renormalized strange quark masses and tnf is the mean of the 
pion masses before and after the reweighting, and V is the physical volume [4]. The uncertainty in 
the volume scaling is currently under investigation. 

The cost for estimating the reweighting factors Wjii), Wjiii), W/i+iv) is roughly independent of 
the ensemble parameters. It is fixed by demanding the stochastic noise to be much smaller than 
the fluctuations [5]. The error bars in Figure 1 indicate the stochastic error for each configuration. 
The total number of Gaussian noise vectors was 48 for W/ii) and W(iii), and 24 for W/i+iv). Since 
the stochastic error is large in Figure 1 these numbers have to be increased in an application of the 
method. 

The numerical cost per configuration (core hours) is ~ 30% of one trajectory of length 2 MDU. 
But since the autocorrelation times of observables are typically much larger than this, reweight- 
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ing is only needed every n = v/2 S> 1 trajectories, effectively cutting the cost to ~ (30/n)%. 

7. Outlook 

In the first numerical test presented in the last section a value for 7 * is used that was determined 
in another setup and at different quark masses. A dedicated tuning needs only a small number of 
configurations and has the potential to lower the fluctuations significantly. Furthermore the cost 
per configuration can be reduced by combining the reweighting factors W(ii), W(iii) into one estima¬ 
tor. Finally, the mrw-package will be upgraded to openQCD-1.4, including support for periodic 
boundary conditions. 

Acknowledgements: This work was funded by the Deutsche Forschungsgemeinschaft (DFG) in 
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